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Abstract 

This paper aims at setting the keystone of a prospective theoretical study 
on the role of non- monotone interactions in biological regulation networks. 
Focusing on discrete models of these networks, namely, Boolean automata 
networks, we propose to analyse the contribution of non- monotony to the 
diversity and complexity in their dynamical behaviours. More precisely, in 
this paper, we start by detailing some motivations, both mathematical and 
biological, for our interest in non-monotony, and we discuss how it may 
account for phenomena that cannot be produced by monotony only. Then, 
to build some understanding in this direction, we propose some preliminary 
results on the dynamical behaviour of some specific non-monotone Boolean 
automata networks called XOR circulant networks. 

Keywords: Discrete dynamical systems. Boolean automata networks, 
non-monotony, dynamical behaviours. 



1. Introduction 

The introduction of Boolean automata networks by McCulloch and Pitts 
in m and Kauffman in [21 [3] has initiated many developments in the study 
of discrete dynamical systems at the frontier of biology, mathematics and 
theoretical computer science. In the context of modelling biological regu- 
lation networks, the pertinence of abstract networks was deeply motivated 
by Hopfield and Kauffman in the respective contexts of neural and genetic 
networks. Among other things, Hopfield showed in [U [5] that threshold 
Boolean automata networks allow to highlight the fundamental neural con- 
cepts of associative memory and learning. In [SI [7j, on the basis of the 
breakthroughs of Jacob and Monod [HI [9] , Kauffman put emphasis on the 
Boolean nature of genes that are simply either actively transcribing or not. 
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These works as well as Thomas' |10^ ITU [T2] placed formal approaches at 
the centre of the understanding of dynamical behaviours and complexity in 
biology. In particular, both these works claimed that theoretical frameworks 
would certainly allow biologists to bypass the observational knowledge which 
cannot, alone, lead to general conclusions. Since then, numerous theoretical 
studies have been carried out to acquire a better understanding of these 
networks, from the computational complexity standpoint [T3 l [T ^ [T5l [T6t [T7] 
as well as from the standpoint of the characterisation of their dynamical 
behaviours [IHlIISlEOlEIlEaESlEllEaES]. 

In the lines of these studies and in order to complete them, we propose 
in this paper to tackle the question of the role of non-monotony in Boolean 
automata networks. This question seems to be missing in classical literature 
dealing with Boolean automata networks as models of biological networks 
and, in particular, as models of genetic regulation networks. Indeed, on 
the one hand, the underlying interaction structure of Boolean models of 
genetic regulation networks are often represented by signed digraphs where 
vertices represent genes and arcs, which are labelled either by a plus or 
a minus sign, represent directed actions of genes on one another, either 
activations or inhibitions. This way, a gene that tends to influence the 
expression of another gene is supposed to be either one of its activators 
or one of its inhibitors. It cannot be both. That is, it cannot act as an 
activator under some circumstances and act as an inhibitor under some 
others. This interpretation of gene regulations leads to define monotone 
Boolean automata networks as studied in [271 [28l [29l l30l [HTj l32] and [331IM1 
[35| [361 157] from theoretical and applied points of view respectively. On the 
other hand, the class of linear networks has also been studied. This class 
contains in particular the special non-monotone networks in which all local 
functions are XOR functions. In [38j, Cull based his study on [U |39l 
and developed an algebraic description of the dynamical behaviour of linear 
networks. In [41j, Snoussi gave a characterisation of behaviours of very 
specific XOR networks. But the global dynamical properties of general non- 
monotone networks have not yet been studied nor has the impact of non- 
monotone interactions yet been examined per se. 

Our recent studies on Boolean automata networks have however brought us 
to believe that non-monotony may be one of the main causes of singular 
behaviours of Boolean automata networks. Thus, this research axis seems 
very pertinent in the context of biological regulation networks. The present 
paper provides the grounds of a prospective study on non-monotony in net- 
works. In this context, we develop two lines. First, with some examples, we 
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give some insights that support the importance of non-monotony and the 
idea that it may be responsible for pecuhar network dynamical behaviours. 
Second, to serve as a tangible starting point and build intuition, we present 
some primary results concerning a particular class of non-monotone net- 
works called XOR circulant networks. 

In Section [2| we provide general definitions and notations about Boolean 
automata networks that are crucial for the sequel. Section [3] details why we 
believe that non-monotony is in some sense at the centre of the existence 
of specific dynamical behaviours. Section |4] presents preliminary results 
concerning XOR circulant networks. In particular, it gives properties of 
their trajectorial behaviours by focusing on convergence times and of their 
asymptotic behaviours by characterising attractors. Eventually, Section [5] 
proposes perspectives to this first work on non-monotony. 

2. Preliminary elements on Boolean automata networks 

Informally, a Boolean automata network involves interacting elements whose 
states, which either equal (inactive) or 1 (active), may change over time 
under the influence of the states of other network elements \18\ H2] . This 
section formalises this description by presenting the main definitions and 
notations which are used in this paper. 

2.1. Structure and local transition functions 

A Boolean automata network N of size n is composed of n elements called 
automata which are, by convention here, numbered from to n — 1. For any 
automaton i oiV = {0, . . . , n — 1}, the set of possible states Xi of i is {0, 1}. 
Let us assume that the time space T is discrete, i.e., T = N. A configuration 
of corresponds to the allocation of a value of {0, 1} to every automaton 
of N . It can thus be represented by a vector x = {xq, . . . ,Xn-i) S {0, 1}" 
and {0, 1}" is then the configuration space of A^. Abusing language, we 
will denote by x{t) (resp. Xi{t)) the configuration of A^ (resp. the state 
of automaton i) at time step t £ T- Given an arbitrary configuration x G 
{0, 1}", the density of x is defined as d{x) = ^ • \{xi \ {i G V) A {xi = 1)}|. 
In our context, we focus particularly on switches of automata states starting 
in a given network configuration. For this reason, the following notations 
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fo(x) = -^xo V {xi A X2) 

fl(x) =XoV X2 
f2{x) = ^Xo 

Figure 1: (left) An interaction graph of a Boolean automata network of size 3 and (right) 
the local transition functions of its automata. 

for network configuratfons will be useful: 
Vx = (xo,...,x„_i) G {0,ir,Vi gV = {0,...,n-l}, 

X — (^^O) • • • ; ^i— 1) • • • ) 2Jn— 1) 

and, yW C V, = W\ (1) 

Thus, in particular, 0* where i £ V (resp. 0^ where W C V) denotes 
the network configuration in which automaton i has state 1 (resp. all the 
automata belonging to W have state 1) and all other automata have state 0. 
The underlying interaction structure of N can be represented by a digraph 
G = {V, A), called the interaction graph of N. In this digraph, V equals the 
set of automata of A'^. A CI V x V is the interaction set. For any automata 
i,j £ V, it satisfies {j,i) G A if and only if j effectively influences i, that is, 
in some network configurations (but not necessarily in all of them), the state 
of j may cause a change of states of i (see Equation[2]below). As an example. 
Figure [1] (left) pictures the interaction graph of a Boolean automata network 
of size 3, where A = {(0, 0), (0, 1), (0, 2), (1, 0), (2, 0), (2, 1)}. Interaction 
graphs specify what influences apply to each automaton of a network A'^. 
The nature of these influences are defined by the local transition functions 
fi : {0, 1}"" — )• {0, 1} which are associated to each automaton i of N (as in 
Figure [1] (right)) such that: 

3x G {0, 1}", Mx) / Mx^) ^ (j, i) G A. (2) 

Thus, a Boolean automata network is entirely defined by the set of local 
transition functions of its automata. 

2.2. Updating modes and transition graphs 

To determine the possible behaviours of a network, it remains to be specified 
how automata states are updated over time. The most general point of 
view consists in considering all possibilities. That is, assimilating networks 
with state transition systems, in each configuration, 2" — 1 transitions are 
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considered, one for each non-empty set of automata whose states can be 
updated. More precisely, yW / C 1/, we define the update function 
Fw : {0, 1}" ^ {0, 1}" such that: 

yx G {0,ir,Vi G V, Fw{x)^ = {^'^''^ if ^ e 

I Xi otherwise. 

Then, according to the most general updating mode, the global network 
behaviour is given by the general transition graph Qg = ({0, l}",Tg) where 
Tg = {{x,Fw{x)) I X e {0,1}", ^ C y} [IHlllllllS]. In this graph 
which usually is a multigraph, arcs can be labelled by the set W of automata 
that are updated in the corresponding transition {x,Fw{x)). For the sake 
of clarity, in the examples of this paper, arcs with identical extremities are 
represented by a unique arc with several labels. Figure [2] (top) depicts the 
general transition graph of the network presented in Figure [T| 

Transitions (x, Fi{x)) that only involve the update of one automaton i £ V 
are called asynchronous transitions. Transitions {x, F\y{x)),\W\ > 1 that 
involve the update of several are called synchronous transitions. The sub- 
graph Qa = ({0,1}", Ta) of Gg whose set of arcs Ta = {{x, F^iy{x)) I X G 
{0,1}", i £ V} equals the set of asynchronous transitions of the network 
is called the asynchronous transition graph. This graph defines the asyn- 
chronous updating mode according to which, in each configuration, only n 
transitions are considered, one for each automaton that can be updated 
alone. This updating mode has been widely used in studies of Thomas and 
his co-workers in jll[ [T^ [2^ l^^l ITTl [53] . An illustration of an asynchronous 
transition graph is given in Figure [2] (bottom). 

Because both the general and the asynchronous transition graphs are very 
large graphs, in some cases, to draw some intuitions, it may be necessary 
to restrict our attention to the transitions that are allowed under a specific 
deterministic updating schedule u. This amounts to considering a transition 
graph Gu = ({0, l}",r„) which is the graph of a function F[u] : {0, 1}" — >• 
{0,1}" (i.e., Tu = {{x,F[u]{x)) \ x G {0,1}"}). This function has the 
following form: F[u] = i^Wp_i ° ■ . . ° Fyi/^ o F^^ where p G N and V/c < 
C V . It is called the global transition function associated to the 
updating schedule u, that updates simultaneously all automata in Wq, then 
updates simultaneously all automata in Wi . . . This second point of view has 
been adopted in [381 HHl [50l [SU [52l [53] following the introduction of hlock- 
sequential updating schedules by Robert in [181 [5l] (see Figurejs]). Section[4]is 
set in similar lines. It focuses on the parallel updating mode tt which consists 
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Figure 2: (top) General and (bottom) asynchronous transition graphs of the Boolean 
automata network of Figure [T] 

in deterministically updating all network automata at once in each network 
configuration. In this case, the global transition function is F[Tr] = Fy such 
that yi E V, F[7r](x)i = fi{x) and the network behaviour is considered 
to be described by the graph of -F[vr], that is, the transition graph Q-^ = 
({0,1}", T,) where = {{x , F[ii]{x)) \ x G {0,1}"}. 

2.3. Dynamical behaviours and non-monotony 

Let be an arbitrary Boolean automata network. Consider any updating 
mode u among those defined above and let Qu be the corresponding tran- 
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(a) (b) (c) 

Figure 3: Transition graphs of the Boolean automata network of Figure [l] associated to 
(a) the parallel updating mode, (b) a block-sequential updating schedule u whose global 
transition function is F[u] — -F{i.2} o F^oy and (c) a particular block-sequential updating 
schedule s, called sequential, whose global transition function is F[u] = F^iy o _F{2} ° ^{o}- 
Recurrent configurations appear in grey. 

sition graph. Let x G {0, 1}" be a configuration of N. We call trajectory 
of X any path in that starts in x. Terminal strongly connected, compo- 
nents of Qu are called the attractors of and constitute the asymptotic 
behaviours of A^. Their size equals the number of configurations that they 
contain. Configurations that belong to an attractor are called recurrent con- 
figurations. Attractors of size 1 (resp. of size strictly greater than 1) are 
called stable configurations (resp. stable oscillations). When G = Gu is the 
transition graph associated to a deterministic updating schedule u, stable 
configurations correspond to fixed points of the global transition function 
F[u] and stable oscillations of size p, which are rather called limit cycles of 
period p in this case, correspond to oriented cycles in Gu- As an example, 
Figure [3] shows that the network of Figure [T] admits one unique attractor, a 
limit cycle of period 2, under any of the three deterministic updating sched- 
ules considered. The precise definition of this limit cycle, however, differs 
in each case. In particular, as proven in \50\ 155). no configurations besides 
stable configurations are recurrent under the parallel updating schedule as 
well as under the sequential updating schedule which updates one automa- 
ton at the time been. Furthermore, Figure [2] shows that the same network 
admits one unique attractor, a stable oscillation of size 8, indifferently when 
it is subjected to the asynchronous or general updating modes. 

By analogy with continuous functions, the local transition function fi of an 
automaton i gV is said to be locally monotone in j G V if, either: 

Vx= (xo,...,x„_i) G {0,1}", 

fi (-^0) ■ • • ) -^j— 1 ) 0) ) ■ ■ • ) -^n— 1 ) ^ fi (3^0 ) • • ■ ) "^J — 1 ) 1) "^ji-l-l ) • • • ) ■^n—l) 
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or: 

Vx = {xo,...,Xn-l) G {0,1}", 
fi (-^O) • • ■ 1 ■^j—lj 0, Xj-\-i ) ^ fi (^Xq, • • • , Xj—lt 1 ) -^J+l ) • ■ ■ ) ''^n— l) • 

In other terms, fi is locally monotone in j if, in the conjunctive normal form 
of fi{x), either only xj appears or only -iXj does. The function fi is said 
to be locally monotone or simply monotone if it is locally monotone in all 
j £ V. It is said to be non (locally) monotone otherwise. In this latter 
case, there exists j £ V such that in some configurations, the state of i 
tends to imitate that of j and in some other configurations, on the contrary, 
the state of i tends to negate that of j. When all functions fi, i G V, are 
monotone, the network is said to be monotone itself. Otherwise, if at least 
one local transition function is non-monotone, the network is said to be 
non-monotone. 

3. Motivations 

To put forward the importance of studying non-monotony in discrete models 
of regulation networks, let us first recall the fundamental concept of genetics 
establishing that a gene is a portion of the dna which is transcribed into a 
mRNA (the gene is then said to be expressed) that is itself translated into one 
or several proteins, called the products of that gene. Because proteins can 
influence the transcription and translation stages, genes have the possibility 
of interacting with one another through their products. Further, because the 
effect of a protein may depend on its concentration in the cell, genes may 
have different effects on one another. If gene gj influences the expression 
of gene gi via one of its protein products p, then it may do so differently 
according to the concentration of p in the cell. As an illustration, let us 
consider the infection of a bacterium Escherichia coli by a phage A [56] at 
high temperature. The genetic regulations that allow a phage A to enter 
its lysogenic and lytic cycles^ involve two genes, Cro and cl. Indeed, at high 
temperature, Cro influences itself and the nature of this influence is different 
according to the concentration of its protein product pcro- In IS3I5H], the 
authors show that if the pcro concentration in the cell is low (resp. high), 
Cro tends to activate (resp. inhibit) its own expression, whatever the Pd 



^The lysogenic cycle of a phage A is the stage where its genome is inserted in the 
genome of the bacterium. Its lytic cycle is the stage where it replicates, leading in fine to 
the death of the bacterium. 
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concentration is. This induces an increase (resp. a decrease) of the pcro 
concentration. Supposing that automata networks are reasonable models of 
genetic networks, such duality in the influence of one gene gj on the state 
of the same or another gene gi corresponds precisely to the formal notion of 
non- monotony, specifically, that of fi{x) with respect to xj. 

In the light of some recent developments j43l |59] , the study of non- monotony 
in discrete networks also seems pertinent from a different, mathematical 
standpoint. Indeed, first, let us consider the problem of modularity in gene 
regulation networks which is essential in the context of biology. Modules 
are informally defined as independent groups of interacting genes. More 
precisely, they involve minimal sets of genes (or biobricks) that own indepen- 
dent behaviours specific to real biological functions. Until now, the notion of 
modularity used in the literature relies on structural parameters [601 EU [62] . 
For instance, modules are often defined simply as strongly connected com- 
ponents. However, although such structural definitions are natural, they 
can only lead to structural results, necessarily failing to reveal biobricks. 
In |59j , a new notion of modularity is introduced for asynchronous multi- 
state automata networks considered as discrete models of gene regulation 
networks. It is supported by dynamical considerations that take into ac- 
count the asymptotic behaviours of these networks. To detail this, let us 
consider a network N whose associated interaction graph is G = (y,A). A 
modular organisation of a network is defined in [SH] by an ordered partition^ 
(Vo, . . . , Vfc_i) of the set of network automata. The subgraph of G induced 
by Vi is said to be a module of N if the ordered sets Vi, i < k, satisfy 
£ = Sq . . . £k-i, where £w, W C V denotes the asymptotic behaviour 
of the subnetwork of induced by W and is a composition operator de- 
fined in |59]. In this context, the authors of |59j show that any topological 
ordering^ on the set of strongly connected components of the underlying 
structure of a network N does indeed define a modular organisation of A^. 
However, in the general case, strongly connected components are not mini- 
mal modules and thus do not allow to reveal biobricks. In some cases, they 
can be decomposed into smaller independent sub-modules. What is inter- 
esting is that for Boolean automata networks, all encountered examples of 



^An ordered partition {So, • • • , Sk-i) of an arbitrary set 5*, is defined hy k < \S\ non- 
empty ordered subsets Si C S such that 5 = 1+]^;, Si {i.e., Wi,j < k, i ^ j ^ SiH Sj — 

^If G = {V, A) is a digraph, a topological ordering of G is a linear ordering of the 
vertices in V such that, V(i, j) £ A, i comes before j in the ordering. 
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0®- KSO ^<1' ^ {x^{xo®xi),x^^{xo®xi)} 



{0},{1} {1} {0},{1},{0,1) 11} 




Figure 4: Top panel: Generic description of the four smallest Boolean automata networks 
that satisfy the conditions of Proposition [l] Bottom panels: (left) generic (see proof of 
Proposition IT]) asynchronous and (right) general transition graphs of these networks. 



non-decomposable strongly connected components involve non-monotony. 

In different lines, the importance of non-monotony can also be seen by 
adapting some results presented in [23]. This produces Proposition [T] be- 
low which relates non-monotony to non-trivial changes in the dynamical 
behaviours of a network when synchronism is added to its asynchronous 
behaviour. In this proposition, a transition {x,F\y{x)) is said to be se- 
quentialisable if there exists a series of consecutive asynchronous transitions 
(x, F^iy{x)), Fjj} o . . . that start in x and end in Fwix). 

Proposition 1. The smallest Boolean automata networks that have non- 
sequentialisable synchronous transitions and significantly different limit be- 
haviours under the asynchronous and general updating modes are non-mono- 
tone. 

Proof. Let us find the smallest network with a non sequentialisable syn- 
chronous transition. Obviously, this network needs to have more than one 
automaton and if it has size 2, then, to have a non-sequentialisable syn- 
chronous transition, its general transition graph needs to contain a subgraph 
of the following form: 

{1} 



{0} 

where x^'^ = x^^'^^ = x^'' (see Equation 1). Moreover, to have significantly 
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different asymptotic behaviours under the asynchronous and general up- 
dating modes, the synchronous transition (x,x^'^) must go out of a set of 
configurations that induces a stable oscillation under the asynchronous up- 
dating mode. Thus, the general transition graph of N must have the form 
of the general transition graph pictured in the bottom right panel of Fig- 
ure |4j Then, only two functions /o are possible. If in configuration x above, 
xq = 1, then, fo{x) : x i— )■ xq © xi where ® denotes the XOR connector^. 
If in configuration x above, xq = 0, then fo{x) : x i— )• -■(xq ® xi). The 
function /i is defined similarly. In conclusion, there are four smallest net- 
works satisfying the properties of Proposition [T| They have size 2 and their 
interactions graph equal the graph pictured in the top panel of Figure |4j 
Their two local interaction functions /o and /i either equal x i— )• xq ® xi or 
X I-)- -i(xo ® Xi). □ 

The reasons evoked in this section that led us to focus on non-monotony 
in automata networks emphasise the apparent importance of non-monotone 
functions. In the next section, to initiate an analysis of the behaviours of 
general non-monotone networks and develop some intuition in this direction, 
we focus on a specific class of non-monotone networks, namely xOR circulant 
networks, and study some of their dynamical properties. 



4. XOR circulant networks 

Before we present some results on the trajectorial and asymptotic dynamical 
behaviours of xOR circulant networks, let us first introduce some definitions 
and preliminary properties in relation to these. 



4-1. Definitions and basic properties 

A circulant matrix C is a matrix of order n whose i^^ row vector Ci {i < n) 
is the right-cyclic permutation with offset i of its first row vector Co so that 
C has the following form: 



/ Co Cl C2 

Cn-l Co Cl 
Cn-2 C„_i Co 



V Cl 



C2 C3 



Cn-l\ 

Cn-2 

Cn-3 

Co / 



"Va, 6 G {0, 1}, a e fe = (a A ^6) V (^a A h). 
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For any integer k > 2, a k-xOR circulant network of size n > A; is a Boolean 
automata network with n automata that can be numbered so that the fol- 
lowing four properties are satisfied: (?) the adjacency matrix C of the net- 
work interaction graph G = {V^A) (i.e., the n x n matrix C defined by 
Vi, j G V, Cij = 1 <^=^ (j, i) G A), called the interaction matrix for short, 
is a circulant matrix, (ii) each row Ci, i & V, of this matrix contains ex- 
actly k non-null coefficients (i.e., Vi G V, X^jgyCjj = degQ{i) = k), (iii) 
Co,n-i = Cn-i = 1 and (iv) the local transition function /j of any automaton 
i G V is a XOR function: 

Vx G {0, 1}", fiix) = 0C,j ■xj = Y^ Ci,j ■ X, [2], 

where, for any integers a and 6, a \b] stands for a (mod b). In the sequel, for 
the sake of simplicity, xOR circulant networks are considered to be subjected 
to the parallel updating mode so that if x = x{t) G {0, 1}" is the network 
configuration at time step t G T, then the network configuration at time 
step t + 1 equals x{t + 1) = F{x) = C ■ x (where operations are supposed to 
be taken modulo 2). Thus, a xOR circulant network is completely defined 
by its interaction graph G = {V, A) or by its interaction matrix C. Figure [s] 
pictures two examples of 2-xOR circulant networks of size 5. Let us note that 
one of the four networks satisfying Proposition [T] and defined in Figure |4] is 
also a 2-xOR circulant network of size 2. 

Let us note that by point (iii) in their definitions, A;-xOR circulant networks 
have Hamiltonian circuits underlying their structures. When automata are 
numbered as suggested in this definition, these circuits are composed of the 
set of arcs {(i, i + 1 [n]) | i G V} C A. More generally, it can be shown that 
each non-null coefficient Cj = Cqj of a circulant interaction matrix C induces 
gcd (n, j) independent circuits of length n/gcd (n, j) in the interaction graph 
G of the corresponding network. Now, it has been shown that to have several 
stable configurations and/or stable oscillations. Boolean automata networks 
need to have circuits underlying their interaction graphs [HI El [63]. Thus, 
the presence of circuits underlying the structures of k-XOR circulant networks 
guarantees that these networks have interesting, non-trivial dynamical be- 
haviours. 

Any k-XOR circulant network N can be seen in terms of cellular automata. 
Indeed, if has size n and interaction graph G = {V, A), it can be modelled 
by the finite one- dimensional cellular automaton that has n cells assimilated 
to the n automata of N and that satisfies what follows. The neighbour- 
hood of a cell i G V equals the in-neighbourhood of automaton i in N: 
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/o(x) 


= xo e 


B Xi 













1\ 


fl(x) 


= Xoi 


Bxi 


1 


1 













= Xli 


dX2 





1 


1 








fsix) 


= X2i 


BX3 








1 


1 





fiix) 


= Xsi 


B Xi 


\o 








1 


1/ 



fo{x) = X2 S 
= Xo S 
f2{x) = Xl < 

fsix) = Xo i 

fi{x) = Xl < 



! X4 
! Xi 
! X4 
! X2 
) X3 



^0 1 1\ 

10 10 
10 1 
10 10 

Vo 1 1 o/ 



Figure 5: Interaction graphs, local transition functions and interaction matrices of two 
2-xOR circulant networks of sizes 5. 

M ={j eV \ e A}. The local rule 7 : {0, ^ {0, 1} of the cel- 
lular automaton is defined similarly to the local transition functions of A^: 
7((x^)^g_V') = ®£(=j\f^e- III the sequel, we use this formalisation to exploit 
tools drawn from the theory of cellular automata. Thus, if x = x{0) G {0, !}"■ 
is an initial configuration of N, we consider the corresponding space-time 
diagram, that is, the grid of {0, 1}" x T whose line t G T represents x{t), i.e., 
the configuration of N at time step t. The trace of cell or automaton i gV 
then corresponds to column i of this grid, that is, to the sequence {xi{t))tQj-. 
Also, for an arbitrary configuration x G {0, l}'^ and an automaton i e V, 
Si{x) denotes the configuration that satisfies Vj G V, Si{x)j = X2i-j \n]- 
It is called the symmetric of x with respect to i. We write N to denote 
the symmetric of N, that is, the k-XOR circulant network whose interaction 
matrix is *C In the sequel, by default, M^(i) (resp. A/'^(i)) denotes the in- 
neighbourhood (resp. the out-neighbourhood) of automaton iinN and A/'~(i) 
(resp. J\f~^{i)) denotes its in-neighbourhood (resp. its out-neighbourhood) 
in A^. This way, for any two automata i,j G y, j G J\f^{i) <^=^ j G J\f^{i). 
The global transition function of N is denoted by F if that of N is denoted 
by F. One last convention that is used throughout the sequel is the fol- 
lowing. By default, unless N is the symmetric of another fc-XOR circulant 
network that was introduced before, its automata are supposed to be num- 
bered as suggested above in the definition of /c-xOR circulant networks so 
that c„_i = Co,n-i = 1- This way, {{i,i [n]) | i G F} C A defines a 
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Hamiltonian circuit in the structure of N and {(i + 1 [n],i) \ i £ V} C A 
defines a Hamiltonian circuit in the structure of its symmetric A^. 

To end this paragraph, we list some basic properties of XOR circulant net- 
works that follow directly from the definitions of XOR functions and circular 
matrices: 

Proposition 2. 

1. The number of k -XOR circulant networks of size n equals (^Z}) • 
Any fc-XOR circulant network of size n satisfies the following properties: 

2. Configuration (0, . . . , 0) is a stable configuration. 

3. Configuration (1, . . . , 1) is an antecedent of (0, ... , 0) if k is even or a 
stable configuration if k is odd. 

4. The trajectory of a configuration x is isomorphic to that of any con- 
figuration y which is a circular permutation of x . 

4.2. Results 

4.2.1. General k-xOR circulant networks 

First, in this paragraph, we concentrate on general A;-XOR circulant networks 
and exploit the cellular automata formalisation presented above to derive 
some features of the dynamical behaviours of these networks. 

Lemma 1. Let N be a k-XOR circulant network of size n with automata set 
V and symmetric global transition function F. For any automaton i £ V, let 
Mi{t), t €zT, denote the set of automata which have state 1 in configuration 
F\W). Then, Vx(0) G {0, 1}", Vt G T, = e,eA/,(t) ^i(O)- 

Proof. We prove Lemma [l] by induction on t G T. 

For t = 0, Mj(0) = {i} holds by definition of configuration 0*. Thus, Vx(0) G 
{0,1}", ^.(0) = e,eAf.(0)^.(0)- 

Now, suppose that Va;(0) G {0, 1}"", Xi{t) = Xj{0) and consider the 

initial configuration y(0) G {0, 1}". 

Since y{t + 1) = F*"'"^(y(0)) = F*(y(l)), applying the induction hypothesis 
to configuration x(0) = y{l) yields: 

j6M,(t) 
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By definition, Vj G V, yj{l) = /,(y(0)) = e^e^^-y) y^O) = e,e;^+(,) ^^0)- 
Thus, with the commutativity and associativity of the operator, we can 
derive that: 

y^{t + l)= ( ym) 

ym 

{e \Af-{i)nM,{t)\=i [2]} 

Now, let us remark that Vt G T, F(0^^'^*^) = g*^''-*^^'' by definition. Then, 
G y, of = 1 if and only if 1X^^(0 n Af,(t)| = 1 [2]. From this follows 
yi{t + 1) = ®jeAU{t+i) Vji^) and then Wt G T, = 0jeAf,(t) ^i(O)- ^ 

Lemma 2. Lei N be a k-xOR circulant network of size n with automata set 
V and global transition function F. For any automaton i £ V , and for any 
configuration x G {0, 1}", it holds that F{Si{x)) = Si{F{x)). 

Proof. For any j G V, the following holds: 

F{Si{x))j = {Si{x))i = X2i-e [n] 



{£ s.t. 2i-e [n] 



{e s.t. jeAf-(2i-e [n])} 

If j G M~{2i — i [n]), then all automata 1,1' G ^ of such that / — 
I' = j — {2i — £) [n] are such that / G M~{1'). In particular, if automaton 
3 G J\f-{2i -l[n]), then I G 7^"(2i - j). Hence: 



{^s.t. jeM~i2i-i [n])} l<^M-{2i-j [n]) 

= F{x)2i-j 

= {smm^ 

and Lemma [2] follows. □ 

Proposition 3. Let N be a k-XOR circulant network of size n with automata 
set V and global transition function F. For any automaton i £ V and for 
the initial configuration x{0) = 0*, it holds that \/t G F, F^{x{0)) = Si{x{t)). 
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Proof. Proposition [3] is proven by induction on t £ T- Let t = 0. Property 
F*(rc(0)) = Si{x{t)) is true because x(0) = 0*. Suppose that it is true 
for t G T. JThen, we have F^+^{x{0)) = F(F*(x(0))) = F{Si{x{t)). By 
Lemma [ij F{Si{x{t)) = Si{F{x{t)) = Si{x{t + 1)), which is the expected 
result. □ 

Remark that this result is due to the fact that F and F are the global 
transition functions of two symmetric k-XOR circulant networks that are 
isomorphic by definition (see Figure [6]). Proposition [s] implies that, for any 
automaton i £ V, the space-time diagram of {0''{t))t<=T is the symmetric 
space-time diagram of (0^'^*'')tg7- with respect to i and is related to the 
trace of automaton i. Thus, the space-time diagrams of configurations of 
density ^ carry information on the global behaviours of N. We examine 
further these properties in the following results. 

Proposition 4. Let N be a k-xOR circulant network of size n with automata 
set V and global transition function F. The maximum convergence time, 
i.e., the maximal transient trajectory length, is reached by configurations 
of density ^. Moreover, let be the period of the attractors reached by 
configurations of density ^ . Then, for any configuration x of N, the period 
of its attractor (i.e., of the attractor that is reached by the network when it 
is initially in configuration x ) divides p-t . 

Proof. Since all configurations of density ^ are cyclic permutations of one 
another, by Proposition [2] 4 they all have isomorphic trajectories so that they 
all hit their limit set at the same time t-t and they all have the same period 
p*. Now, consider configuration x and automaton i. By Proposition [3j the 
space-time diagram of (0*^'*^*'')tg7- is the symmetric space-time diagram of 
((f (t))(g7- with respect to i. Thus, the space-time diagram of (0*^'^*^)(g7- 
hits its limit set at time and its period is p-t. This means that, "ii £ N , 
the trace of automaton i has period p^ and hits its limit before t*. Thus, the 
trajectory of x reaches its limit set before t* and its period divides p*. □ 

4.2.2. 2-XOR circulant networks 

Let us now concentrate on 2-xOR circulant networks of arbitrary size n 
and pay particular attention to the space-time diagrams of configurations 
of density ^. We define the interaction- step of such a network N as the 
smallest integer s ^ 1 < n such that Mi EV,{i,i + s [n]) G A. As illustrated 
in Figure [g] (a) and (6), when s = the space-time diagram is the Sierpinski 
triangle. For other values of s, space-time diagrams seem like deformed 
Sierpinski triangles. From these observations results the following lemma. 
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Lemma 3. Let N be a 2-xOR circulant network of size n with interaction- 
step s = 0. The following holds: 

yi evyqe N, x,(2«) = x(i„2'^) n(0) ® 2;i(0). 
Proof. Lemma [3] is proven by induction on q. 

Let i G y be an arbitrary automaton and let q equal 1 initially. Then, 
obviously, the following is true: 

Xi{2) = X(^i_i) [„](!) 0x^(1) 

= a^(i-2) [n](0) ex(i_i) [„](0)exi(0) 

= X{i~2) [n]{0) ®Xi{0), 

and the basis of the induction holds. 

Now, let us assume as induction hypothesis that Xi{2'^) = a;(j_29) [„](0)©Xj(0) 
is true for q £ N. In the sequel, we pay particular attention to states 

a = Xi{0), b = a;(j„2'j-i) [n](0), c = X(^i_2q) [„](0), d = Xi{2'^~^), 
e = X(j_2<7-i) [„](2«-i) and / = Xi(2«), 

as illustrated in Figure [7j 
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Figure 7: Space-time diagram of an arbitrary 2-XOR circulant network of size n and 
interaction-step s = 0. 



Then, for g + 1, according to the induction hypothesis, we have: 

d = a®b, e = h®c and f = d®e. 
Then, we derive that: 

/ = d®e = (a e 6) e (6 e c) = a e c. 
As a result, we can write: 

Vi G y, Vg G M, Xi{2'i) = x,(0) X(i_2,) [„](0), 

which is the expected result. □ 

We will use this lemma to analyse 2-XOR circulant networks of size n = 2^, 
p G N*, and interaction step s = 0. 

4-2.3. 2-XOR circulant networks of sizes powers of 2 

In this paragraph, we focus on 2-xOR circulant networks of sizes n = 2^, 
where p G N*. Let x = {xo, . . . , Xn-i) G {0, 1}" be a configuration of such 
a network N. We can see x as the concatenation of two vectors of sizes | 
such that x = {x',x"), where x' = {xq, . . . and x" = . . . ,Xn-i) 

are called semi-configurations of x. The repetition degree Sr{x) of x is then 
defined as: 



Sr{x = {x',x")) 



if x' / x", 

5 if {x' = x") A {dr{x') = 6-1) 



Notice that if x = {x',x'), x is said to be a repeated configuration and 
that, in the worst case (i.e., when the repetition degree 6r{x) = log2(n) is 
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maximal), the time complexity of the computation of the repetition degree 
of any configuration x equals n. 

Proposition [5] below characterises the dynamical behaviours of repeated con- 
figurations X G {0, l}'^ of repetition degree 5r{x) > log2(n) — 1 in arbitrary 
2-XOR circulant networks of size n = 2^', p G N*. 

Proposition 5. Let N he a 2-xOR circulant network of size n = 2^, p £ 
N* , and interaction-step s. Configurations x € {0,1}" of repetition degree 
6r{x) > log2(n) — 1 converge towards (0, . . . , 0) in no more than 2 time steps. 

Proof. First, notice that because is a 2-xOR circulant network of size 
n = 2^, p £ N*, there exist only 4 repeated configurations of degree no 
smaller than log2(n) — 1, that is, (0, 1, ... ,0, 1), its dual (1,0, ... , 1,0) and 
(1, . . . , 1) and its dual (0, . . . , 0). Let us consider the two distinct parities 
of s independently. Also, let t G T and let x(t) be either (0, 1, ... ,0, 1) or 
(1,0,. ..,1,0). 

1. If s is even, then, by hypothesis on x{t): 

yi G V, [n]{t + 1) = Xi{t) e X(^i+s-l) [n]{t) = 1- 

2. If s is odd, then, by hypothesis on x(t): 

Vi G V, [„](t + 1) = Xi{t) e [„](t) = 0. 

This, together with Propositions [2] 2 and[2}3, yields the expected result. □ 

From now on, we restrict the study to 2-xOR circulant networks of sizes 
n = 2P, p £ N*, and interaction-steps s = 0. We show that such networks 
necessarily converge towards configuration (0, . . . , 0) in no more than n time 
steps and that initial configurations with an odd number of 1 converge in 
exactly n steps. 

Theorem 1. Let N be a 2-xOR circulant network of size n = 2^, p G N* , and 
interaction-step 0. Any configuration x converges to the stable configuration 
(0, . . . , 0) in no more than n time steps. 

Proof. Since n = 2^, by Lemma [3j we directly draw: 

Vi G V, Xi{n) = Xi{0) e Xi+n [„](0) = Xi{0) Xi{0) = 0. 

This allows to conclude that any configuration x converges to the stable 
configuration (0, . . . , 0) in no more than n time steps. □ 
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Now, let us consider the configurations for which the convergence time is 
maximal. 

Lemma 4. Let N and N' be two 2-XOR circulant networks of respective 
sizes n = 2^^^ and n' = 2^, p £ N* , and interaction-steps 0. Let x' be a 
configuration of size 2^ and x = {x',x') be a repeated configuration of size 
2P+^ Then, for any teT, x{t) = {x' (t) , x' (t)) . 

Proof. Considering an arbitrary repeated configuration x of N, by induction 
on t, we show that Vt G T, x{t) = {x' (t) , x' {t)) . Let us denote by G' = 
(y' , A') the interaction graph of N' . 

By hypothesis, the proposition is true for t = 0. 

Now, consider that x{t) = {x' (t) , x' {t)) for t gT and that 

Vi G V, Xi{t + 1) = e Xi{t). 

Since x{t) is a repeated configuration, we have: 

Vi G V, Xi{t + 1) = © Xi{t) 

= X(i-1+2P) [n](i) ffi a;(j+2P) [„](t) 
= ^(i+2P) [n]{t + !)• 

Thus, x{t + 1) is also repeated and it satisfies: 

Vi G V\ Xi{t + 1) = [.„/] (t) © 

= 4-1) [n'](*)®^UO 

= x',(t + l). 

As a resuh, it holds that x{t + 1) = {x'{t + + 1)). □ 

Proposition 6. Let N he a 2-xOR circulant network of size n = 2^, p G N*, 
and interaction- step 0. Any configuration x such that n ■ d{x) = 1 [n] (with 
an odd number of Is) converges in n time steps exactly. 

Proof. Proposition [6] is proven by induction on p. 

li p = 1, according to Propositions [2j 3 and [5| configurations of repetition 
degree log2(n) — 1 are proven to converge in 2 time steps. Thus, Proposition [g] 
holds for p = 1. 

Suppose that for p = q, any configuration x such that 2'^ • d{x) = 1 [n] 
converges in 2'^ time steps. 

Now, suppose that p = q -\- 1 and consider a 2-XOR circulant network of 
size n = 2'^'^^ and interaction-step 0. Let x be a configuration of size 2"^+^ 
such that n ■ d{x) = 1 [ri\. We show that after 2^ time steps: 
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1. x{2'^) is a repeated configuration of the form x{2'^) = (x'(2'^), x'(2'')). 
By Lemma g Vi G {0,...,2'? - 1}, 2:^(2'?) = Xi(0) X(,+2«) [„](0). 
Hence: 

ViG {0,...,2'?-l}, 

Xi(2«) = X(j+2'J+i) [n](0) ea;(i+2'/) ^(O) = X(^i_^2i) [n](2^)- 

2. x' has an odd number of Is. By Lemma [s] and the property above, 
since Vi G {0, . . .2^ - l},x'i(29) = Xi{2i) = Xi{0) + ^(^+2^) [„](0), each 
automaton of x(0) influences exactly one automaton of x' . If x^(2'') = 
0, then the states of both the automata of x(0) that influence x^(2'') 
must have the same parity. If x^(2'') = 1 then the states of both the 
automata of x(0) that influence x^(2'^) must have different parities. 
Since there is an odd number of Is in x(0), there is an odd number of 
Is in x'(25). 

By Lemma[4| x(2'') behaves exactly hke x'(2''). By the induction hypothesis, 
x' converges in exactly 2' time steps. Thus x converges in exactly n = 2'^"'"^ 
time steps. □ 

5. Conclusion and perspectives 

With this study, we have endeavoured to show that non-monotony is an 
interesting concept per se, despite the lack of speciflc attention it has re- 
ceived so far. On the one hand, to serve as a stepping-stone and acquire 
some initial intuitions in this domain, we have considered a special family 
of non-monotone Boolean automata networks that we named XOR circulant 
networks. In particular, we have focused on the trajectorial and asymp- 
totic behaviours of these networks, considered their convergence times and 
characterised their attr actors. Globally, this preliminary formal analysis re- 
vealed that simple non-monotone networks can exhibit non-trivial, engaging 
properties. On the other hand, more generally and informally, we also have 
put forward several arguments to support the idea that work needs to be 
done to build a better understanding of the role of non-monotony in the 
behaviours of automata networks. In these lines, we have mentioned that 
studies in this context could flnd concrete and relevant applications in bi- 
ology and in particular in the modelling of genetic regulation networks by 
automata networks. In addition, we have given two theoretical arguments 
in favour of our insights by which non-monotony could be responsible for 
singular network behaviours. First, exploiting [59], we have argued that 
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non-monotony may be responsible for the strongly connected components 
of networks being non-separable, minimal functional modules. Second, in 
Proposition [T| with the state transition systems formalism, we have consid- 
ered "synchronism sensitivity" , that is, the property of Boolean automata 
networks to display significant behavioural changes when synchronism is 
added to their automata state updates. And in this context, we have shown 
that the smallest synchronism sensitive Boolean automata networks are also 
non-monotone. 

The issues presented in this paper open many research directions that could 
help develop a better understanding of the precise role of non-monotony in 
formal automata networks and, a fortiori, in real biological regulation net- 
works. One of these perspectives consists in identifying the relations that 
exist between monotone and non- monotone Boolean automata networks. 
In |43j . some preliminary results are derived on synchronism sensitivity. In 
particular it is shown that this property requires specific circuits underlying 
the networks structures. Also, monotone examples of synchronism sensitive 
networks are given. The interesting point is that all of them seem to involve 
a monotone coding of non-monotony. With Proposition [T| this naturally 
raises the question of whether non-monotony (taken in a more general sense 
than what we did formally above) can account in a certain way for the 
synchronism sensitivity in arbitrary monotone and non-monotone networks. 
Thus, Proposition [l] together with the work presented in [43j call for further 
researches in this direction. With sufficient knowledge in this context, we 
then hope to move on to the subject of modularity as developed in [59] and 
work on establishing the exact non-separability conditions of strongly con- 
nected networks. In this context, the first important questions that need to 
be addressed are: "Does there exist monotone strongly connected networks 
that are separable into functional modules?" and "How does non-monotony 
relate to the non-separability of non-monotone networks?" . The relevance of 
these questions lies in that their answers will help understand modularity in 
biological regulation networks, which is a central issue in present biological 
research frameworks such as synthetic biology. Eventually, further analyses 
also need to be done on the dynamical behaviours of XOR circulant networks. 
Indeed, we believe that these networks constitute very promising instances 
of non- monotone networks because of their apparent simplicity and because, 
since they involve underlying structural circuits, their dynamical behaviours 
are potentially diverse and complex. Thus, pursuing in this direction, we 
hope to obtain generalisations of the results that figure above concerning 
the parallel updating mode by relaxing structural constraints step by step. 
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Also, another interesting perspective in this framework is to consider XOR 
circulant networks as state transition systems, under the asynchronous and 
general updating modes. This perspective is motivated in particular by the 
fact that, according to Proposition [T| the smallest synchronism sensitive 
networks are either XOR circulant networks of size 2 and interaction-step 0, 
or networks that have the same structures as these and comparable non- 
monotone interactions. 

6. Acknowledgements 

We thank Florian Rabin for his relevant comments and are indebted to 
the Agence nationale de la recherche and the Reseau national des systemes 
complexes that have respectively supported this work through the projects 
Synbiotic (anr 2010 blan 0307 01) and Meteding (rnsc ai10/11-l03908). 

References 

[1] W. S. McCulloch, W. Pitts, A logical calculus of the ideas immanent in 
nervous activity, Journal of Mathematical Biology 5 (1943) 115-133. 

[2] S. A. Kauffman, Metabolic stability and epigenesis in randomly con- 
structed genetic nets. Journal of Theoretical Biology 22 (1969) 437-467. 

[3] S. A. Kauffman, Homeostasis and differentiation in random genetic con- 
trol networks. Nature 224 (1969) 177-178. 

[4] J. J. Hopfield, Neural networks and physical systems with emergent 
collective computational abilities, Proceedings of the National Academy 
of Sciences of the USA 79 (1982) 2554-2558. 

[5] J. J. Hopfield, Neurons with graded response have collective compu- 
tational properties like those of two-state neurons, Proceedings of the 
National Academy of Sciences of the USA 81 (1984) 3088-3092. 

[6] S. A. Kauffman, Current topics in developmental biology. Vol. 6, El- 
sevier, 1971, Ch. Gene regulation networks: a theory for their global 
structures and behaviors, pp. 145-181. 

[7] S. A. Kauffman, The origins of order, Oxford University Press, 1993. 

[8] F. Jacob, J. Monod, Genetic regulatory mechanisms in the synthesis of 
proteins. Journal of Molecular Biology 3 (1961) 318-356. 



23 



[9] F. Jacob, J. Monod, On the regulation of gene activity, Cold Spring 
Harbor Symposia on Quantitative Biology 26 (1961) 193-211. 

[10] R. Thomas, Boolean formalisation of genetic control circuits. Journal 
of Theoretical Biology 42 (1973) 563-585. 

[11] R. Thomas, On the relation between the logical structure of systems and 
their ability to generate multiple steady states or sustained oscillations, 
in: Numerical methods in the study of critical phenomena, Vol. 9 of 
Springer Series in Synergetics, Springer, 1981, pp. 180-193. 

[12] R. Thomas, Regulatory networks seen as asynchronous automata: a 
logical description. Journal of Theoretical Biology 153 (1991) 1-23. 

[13] P. Florccn, P. Orponen, On the computational complexity of analyzing 
Hopfield nets. Complex Systems 3 (1989) 577-587. 

[14] M. Cosnard, P. Koiran, H. Paugam-Moisy, Complexity issues in neural 
network computations, in: Proceedings of LATIN, Vol. 583 of Lecture 
Notes in Computer Science, Springer, 1992, pp. 530-544. 

[15] P. Koiran, Puissance de calcul des reseaux de neurones artificiels, Ph.D. 
thesis, Ecole normale superieure de Lyon (1993). 

[16] P. Orponen, Computing with truly asynchronous threshold logic net- 
works. Theoretical Computer Science 174 (1997) 123-136. 

[17] A. Gajardo, A. Moreira, E. Coles, Complexity of Langton's ant. Dis- 
crete Applied Mathematics 117 (2002) 41-50. 

[18] F. Robert, Discrete iterations: a metric study. Vol. 6 of Springer Series 
in Computational Mathematics, Springer, 1986. 

[19] E. Coles, J. Olivos, Comportement periodique des fonctions a seuil 
binaires et applications. Discrete Applied Mathematics 3 (1981) 93- 
105. 

[20] E. Coles, Neural and automata networks: dynamical behaviour and 
applications. Mathematics and Its applications, Kluwer Academic Pub- 
lishers, 1990. 

[21] J. Aracena, J. Demongeot, E. Coles, Fixed points and maximal inde- 
pendent sets in AND-OR networks. Discrete Applied Mathematics 138 
(2004) 277-288. 



24 



[22] J. Aracena, J. Demongeot, E. Goles, On limit cycles of monotone func- 
tions with symmetric connection graph, Theoretical Computer Science 
322 (2004) 237-244. 

[23] E. Remy, P. Ruct, D. Thicffry, Graphic requirements for multistability 
and attractive cycles in a Boolean dynamical framework, Advances in 
Applied Mathematics 41 (2008) 335-350. 

[24] J. Demongeot, E. Goles, M. Morvan, M. Noual, S. Sene, Attraction 
basins as gauges of robustness against boundary conditions in biological 
complex systems, PLoS One 5 (2010) el 1793. 

[25] J. Demongeot, M. Noual, S. Sene, Combinatorics of Boolean automata 
circuits dynamics. Discrete Applied Mathematics, in press. 

[26] A. Richard, Local negative circuits and fixed points in non-expansive 
Boolean networks. Discrete Applied Mathematics 159 (2011) 1085- 
1093. 

[27] E. Golcs-Chacc, F. Fogelman-Soulie, D. Pcllcgrin, Decreasing energy 
functions as a tool for studying threshold networks. Discrete Applied 
Mathematics 12 (1985) 261-277. 

[28] M. Cosnard, E. Goles, Discrete state neural networks and energies. 
Neural Networks 10 (1997) 327-334. 

[29] E. Remy, B. Mosse, C. Chaouiya, D. Thieffry, A description of dynam- 
ical graphs associated to elementary regulatory circuits, Bioinformatics 
19 (2003) iil72-iil78. 

[30] C. Chaouiya, E. Remy, P. Ruet, D. Thieffry, Qualitative modelling of 

genetic networks: from logical regulatory graphs to standard Petri nets, 
in: Proceedings of Petri Nets, Vol. 3099 of Lecture Notes in Computer 
Science, Springer, 2004, pp. 137-156. 

[31] O. Colon-Reyes, R. Laubenbacher, B. Pareigis, Boolean monomial dy- 
namical systems. Annals of Combinatorics 8 (2005) 425-439. 

[32] A. S. Jarrah, R. Laubenbacher, A. Veliz-Cuba, The dynamics of con- 
junctive and disjunctive Boolean network models. Bulletin of Mathe- 
matical Biology 72 (2010) 1425-1447. 

[33] L. Mendoza, E. Alvarez-Buylla, Dynamics of the genetic regulatory net- 
work for Arabidopsis thaliana flower morphogenesis. Journal of Theo- 
retical Biology 193 (2) (1998) 307-319. 



25 



[34] L. Mendoza, D. ThiefFry, E. Alvarez-Buylla, Genetic control of flower 
morphogenesis in Arabidopsis thaliana: a logical analysis, Bioinformat- 
ics 15 (1999) 593-606. 

[35] J. Araccna, M. Gonzalez, A. Zuhiga, M. A. Mendez, V. Cambiazo, 
Regulatory network for cell shape changes during drosophila ventral 
furrow formation. Journal of Theoretical Biology 239 (2006) 49-62. 

[36] C. Georgescu, W. J. R. Longabaugh, D. D. Scripture-Adams, E. S. 
David-Fung, M. A. Yui, M. A. Zarnegar, H. Bolouri, E. V. Rothenberg, 
A gene regulatory network armature for T lymphocyte specification, 
Proceedings of the National Academy of Sciences of the USA 105 (2008) 
20100-20105. 

[37] L. Mendoza, F. Pardon, A robust model to describe the differentiation 
of T-helper cells, Theory in Biosciences 129 (2010) 283-293. 

[38] P. Cull, Linear analysis of switching nets. Biological Cybernetics 8 
(1971) 31-39. 

[39] D. A. Huffman, Information theory, Academic Press, 1956, Ch. The 
synthesis of linear sequential coding networks. 

[40] B. Elspas, The Theory of Autonomous Linear Sequential Networks, IRE 
Transactions on Circuit Theory 6 (1959) 45-60. 

[41] E. H. Snoussi, Structure et comportement iteratif dc certains modeles 
discrets, Ph.D. thesis, Universite Grenoble 1 - Joseph Fourier (1980). 

[42] C. Choffrut (Ed.), Automata networks. Vol. 316 of Lecture Notes in 
Computer Science, Springer, 1988. 

[43] M. Noual, Synchronism vs asynchronism in Boolean networks. Tech. 
rep., Ecole normale superieure de Lyon, arXiv: 1104.4039 (2011). 

[44] M. Noual, General transition graphs and Boolean automata circuits, 
submitted. 

[45] M. Noual, S. Scnc, Towards a theory of modelling with Boolean au- 
tomata networks - I. Theorisation and observations. Tech. rep., Ecole 
normale superieure de Lyon and Universite d'Evry - Val d'Essonne, 
arXiv:1111.2077 (2011). 



26 



[46] A. Richard, J.-P. Comet, G. Bernot, R. Thomas' modeUng of biological 
regulatory networks: introduction of singular states in the qualitative 
dynamics, Fundamenta Informaticae 65 (2004) 373-392. 

[47] A. Richard, J.-P. Comet, Necessary conditions for multistationarity in 
discrete dynamical systems. Discrete Applied Mathematics 155 (2007) 
2403-2413. 

[48] J. Demongeot, A. Elena, S. Sene, Robustness in regulatory networks: 
a multi-disciplinary approach. Acta Biotheoretica 56 (2008) 27-49. 

[49] A. Elena, J. Demongeot, Interaction motifs in regulatory networks and 
structural robustness, in: Proceedings of CISIS, IEEE, 2008, pp. 682- 
686. 

[50] J. Aracena, E. Goles, A. Moreira, L. Salinas, On the robustness of 
update schedules in Boolean networks, Biosystems 97 (2009) 1-8. 

[51] J. Aracena, E. Fanchon, M. Montalva, M. Noual, Combinatorics on 
update digraphs in Boolean networks. Discrete Applied Mathematics 
159 (2010) 401-409. 

[52] E. Goles, M. Noual, Block-sequential update schedules and Boolean 
automata circuits, in: Proceedings of Automata 2010, DMTCS, 2010, 
pp. 41-50. 

[53] E. Goles, M. Noual, Disjunctive networks and update schedules, Ad- 
vances in Applied Mathematics, in press. 

[54] F. Robert, Les systemes dynamiques discrets, Vol. 19 of Mathematiques 
k, Applications, Springer, 1995. 

[55] E. Goles, L. Salinas, Comparison between parallel and serial dynamics 
of Boolean networks. Theoretical Computer Science 396 (2008) 247- 
253. 

[56] E. M. Lederberg, Lysogenicity in Escherichia coli strain K-12, Microbial 
Genetics Bulletin 1 (1950) 5-9. 

[57] H. Eisen, P. Brachet, L. Pereira da Silva, F. Jacob, Regulation of repres- 
sor expression in A, Proceedings of the National Academy of Sciences 
of the USA 66 (1970) 855-862. 



27 



[58] D. ThiefFry, R. Thomas, Dynamical behaviour of biological regulatory 
networks - II. Immunity control in bacteriophage lambda, Bulletin of 
Mathematical Biology 57 (1995) 277-297. 

[59] F. Delaplacc, H. Klaudel, T. Melliti, S. Sene, Analysis of modular or- 
ganisation of interaction networks based on asymptotic dynamics, sub- 
mitted. 

[60] R. Milo, S. Shen-Orr, S. Itzkovitz, et al., Network motifs: simple build- 
ing blocks of complex networks, Science 298 (2002) 824-827. 

[61] A. W. Rives, T. Galitski, Modular organization of cellular networks. 
Proceedings of the National Academy of Sciences of the USA 100 (2003) 
1128-1133. 

[62] J. Gagneur, R. Krause, T. Bouwmeester, et al.. Modular decomposition 
of protein-protein interaction networks. Genome Biology 5 (2004) R57. 

[63] A. Richard, Negative circuits and sustained oscillations in asynchronous 
automata networks. Advances in Applied Mathematics 44 (2010) 378- 
392. 



28 



